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Abstract 

The feasibility and design of an acoustic neutrino detection array in the South Pole ice depend on the acoustic 
properties of the ice. The South Pole Acoustic Test Setup (SPATS) was built to evaluate the acoustic characteristics 
of the ice in the 1 to 100 kHz frequency range. The vertical sound speed profile relates to the level of refraction of 
the surface noise and determines the reconstruction precision of the neutrino direction. The SPATS speed of sound 
analysis for pressure and shear waves is presented. 
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1. Introduction 

1.1. Motivation 

A hybrid optical/radio/acoustic neutrino detec- 
tor at the South Pole is described and simulated 
in [J, [2| • The design of an acoustic neutrino detec- 
tion array as part of this hybrid detector requires 
knowledge of the acoustic properties of the detec- 
tor medium. The South Pole Acoustic Test Setup 
(SPATS) 0] was built to evaluate the acoustic at- 
tenuation length, background noise level, transient 
rates and sound speed in the South Pole ice for the 1 
to 100 kHz region. When an ultra-high-energy neu- 
trino interacts in the ice, the hadronic cascade lo- 
cally and instantaneously warms up a cylindrical 
portion of the medium followed by a slow dissipa- 
tion of the heat. The resulting acoustic signature is 
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a pressure wave with a disk-like emission pattern 
perpendicular to the shower direction. The speed at 
which this pressure wave then propagates through 
the South Pole ice is an important factor in event 
reconstruction and transient background rejection. 
In addition to obtain knowledge about the absolute 
value of the sound speed in the ice, it is important to 
map the vertical sound speed profile. The acoustic 
waves are bent towards regions of lower propagation 
speed and the profile dictates the refraction index 
and the resulting radius of curvature. The emission 
disk is deformed more for larger sound speed gradi- 
ents (lower radii of curvature) , making the direction 
reconstruction for the hadronic shower more diffi- 
cult. 

The sound speed depends on the density and tem- 
perature of the medium. From the surface to the 
bulk-ice, in the firn region, the ice becomes gradu- 
ally stiffer due to the increased density. Therefore 
the sound speed rapidly increases in this ice layer. A 
sound speed profile for the firn layer was experimen- 
tally obtained by J.G. Weihaupt 0] in 1963 using 
seismic surface measurements (see also Figure [3]). 
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For larger depths, the density is stable and 
the sound speed can be modeled using a certain 
temperature-dependent ice coefficient as input. The 
temperature of the ice at South Pole has been ex- 
perimentally determined to go from -51°C just be- 
low the surface to an expected -9°C at 2800 meters 
depth, close to the bedrock Therefore, by com- 
bining this temperature profile with the sound speed 
temperature dependence, we can expect a slight 
decrease of the sound speed towards the deeger ice. 
For the in-situ measurements presented in [6], this 
temperature dependence of the speed of pressure 
waves was found to be -2.30 ± 0.17 In l7] a 
certain model for the temperature-dependent ice 
coefficient and J.G. Weihaupt's firn measurements 
are combined to model the pressure and shear wave 
speeds for the South Pole ice. In this model, the 
surface noise is predicted to be mostly trapped in 
the refractive firn layer. The human activity on the 
surface is a likely source of background noise, but 
the sound should not propagate to the region below 
firn. A much smaller and downward refraction for 
the noise originating from the bedrock is expected. 

1.2. Setup 

The South Pole Acoustic Test Setup consists of 4 
acoustic strings. An array of three vertical strings 
(A, B and C) was deployed in the upper 400 m of 
the Antarctic ice cap in January 2007. It was ex- 
tended in the 2007/2008 polar season by the addi- 
tion of a fourth string (D) down to 500 m depth. 
Each string was deployed in the upper part of an Ice- 
Cube [Sj hole and has 7 stages with each one trans- 
mitter and one sensor module. Strings A, B and C 
have stages at 80, 100, 140, 190, 250, 320 and 400 m 
depth while string D has stages at 140, 190, 250, 320, 
400, 430 and 500 m depth. This defines 9 distinct 
instrumented levels. For an overview of the current 
status of SPATS, see 0]. A retrievable transmitter 
or "pinger" was deployed in 6 water-filled IceCube 
holes. The pingeiUJ offers a broadband omnidirec- 
tional emission and has a transmission power of 149 
dB per (/itPa/V) at 1 m distance. It was lowered 
down to a maximum depth of 500 m while contin- 
uously transmitting at 1 Hz. At each SPATS sen- 
sor level, the lowering was stopped for 5 minutes al- 
lowing all sensors to record at least 1 waveform of 
9 seconds at 200 kHz sampling rate. The analysis 
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Fig. 1. Top view of the IceCube detector as of February 
2008, the SPATS (47, 72, 76, 78) and pinger (55, 69, 70, 
71, 76, 77) holes are labeled. Nominal distance between the 
nearest-neighbor IceCube holes is 125 m. 

presented here uses data from geometries where the 
pinger and sensor are at the same depth and 125 m 
apart. Figure 1 shows a schematic of the SPATS 
and pinger geometry. The pinger stage is lowered on 
a ~2700 m long armored cable through which the 
power and trigger are routed to the pressure vessel 
that contains the high voltage (HV) board. Upon re- 
ceiving the trigger signal, a short (~60 //s) HV pulse 
is sent to the piezo-ceramic transmitter that is sus- 
pended about 2 m below the housing. A broadband 
acoustic pulse is then emitted. Both the SPATS ar- 
ray and the pinger are synchronized by GPS so that 
the exact transit times of the pinger pulses can be 
determined. Data were collected with all 9 SPATS 
instrumented levels. An example of a pinger pulse 
recorded by a SPATS sensor is shown in Figure [21 

2. Sound speed analysis 

2.1. Method 

The sound speed is expected to be independent 
of direction in the depths where SPATS modules 
are deployed. In single-crystal ice the sound speed 
depends on the direction of propagation relative to 
the c-axis. But the South Pole ice is polycrystalline, 
with crystal sizes of order 0. 1 cm and with a random 
distribution of c-axis orientations, such that this ef- 
fect averages out [91 . The frequency dependence of 
the sound speed is also negligible. 

The SPATS interstring (transmitting with in-ice 
transmitters and recording on other strings) and 
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Fig. 2. A SPATS sensor waveform and arrival-time (full ver- 
tical line) determination using a fixed threshold (full hori- 
zontal lines). 

pinger data both show shear waves arriving, as ex- 
pected, at roughly double the pressure wave transit 
time for baselines under 200 m. In the case of the 
SPATS transmitters we can understand the shear 
waves as being created at the transmitter-ice bound- 
ary. As for the pinger emitter, it is located in a col- 
umn of water where no shear waves can propagate. 
The observed shear waves are therefore expected to 
be generated at the water-ice boundary where mode 
conversion can occur. In this case the relative shear 
and pressure wave amplitudes depend on the inci- 
dent angle and therefore on the position of the trans- 
mitter in the hole. Arrival times can be extracted 
from the data for both pressure and shear waves and 
for all SPATS instrumented levels. 

The sound speed for each level is defined as the 
distance between pinger and SPATS sensor divided 
by the transit time (time between emission and de- 
tection of the pulse) . The pinger emits a pulse at ev- 
ery start of a GPS second (1 Hz repetition rate). The 
SPATS strings are also synchronized with a GPS 
IRIG-B 100 pps signal so that the start of every sec- 
ond can be determined in the SPATS sensor wave- 
forms. The arrival time (ta) is then defined as the 
time elapsed from the start of the second to the first 
rising edge of the received pinger waveform. This 
edge corresponds to the first sample above a thresh- 
old that depends on the background noise level of 
that specific sensor-channel (Figure [5]) . The delay 
due to cable and electronics (tm) is defined at the 
pinger-sidc as the delay from the start of the GPS 
second to the actual acoustic emission. This delay 



needs to be subtracted from the arrival time in order 
to obtain the transit time (tt) and was measured in 
a laboratory setup to be 1.9±0.1 ms. From this we 
find the transit time: tt=ta-tm- 

2.2. Results 

Three sources of measurement uncertainties are 
taken into account (see Table 1). First, for the 
125 m baseline wc find that the error on t„j comes 
to 0.31% for the 0.1 ms uncertainty. Second, taking 
a ±0.05 ms window for the arrival time assures a 
correct selection of the first peak for the waveform; 
this yields 0.15%. 



Measurement 


Result 


Error on tt 


emission delay 


1.9±0.1ms 


0.31% 


arrival time 


34.00±0.05 ms 


0.15% 


distance 


125.0±0.7m 


0.57% 


Total 




0.66% 



Table 1 

The measurement uncertainties for a transit time of ~32 ms. 

Finally, the uncertainty on the horizontal distance 
between pinger and SPATS sensor needs to take into 
account the uncertainty on the position of the cen- 
ter of each hole. Although the nominal horizontal 
distance between the centers of two adjacent holes 
is 125 m, deviations can indeed be expected due to 
imprecision of initial hole positioning or intentional 
displacement to avoid surface cabling or debris in 
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Fig. 3. VpresBure and Vghear versus depth from SPATS-pinger 
speed of sound analysis, la error bars are indicated. The 
results from J.G. Weihaupt Q, for which no uncertainties 
were cited, are included for comparison. 
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the firn. The position of the center of each hole is 
determined after deployment by a GPS survey with 
a precision of 0.5 m. From there we assume ±0.5 m 
in any direction and 0.57% (0.5\/2 m for 125 m) as a 
conservative error on the distance. After combining 
the three measurement uncertainties, we find the to- 
tal uncertainty to be 0.66% for the 125 m baseline. 
The sound speed analysis can benefit from longer 
baselines since the error on the distance between 
transmitter and sensor is dominant. The minimal 
uncertainty achievable with the 2007-2008 SPATS 
pinger data is 0.1% by going to the maximum base- 
line of 543 m. 

Figure 3 shows all SPATS sound speed data points 
for depths from 80 m to 500 m. There is agreement 
with [j] in the firn region for the pressure wave mea- 
surement. Figure 4 shows a close-up for the 250 to 
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waves and it is consistent with no refraction. Un- 
der the hypothesis of no sound speed gradient below 
250 m depth, the pressure wave speed was measured 
to be 3850 ± 50 m/s and the shear wave speed to be 
1950 ± 50 m/s in the 250 to 500 m depth region. 

3. Conclusion 

We presented the sound speed analysis and results 
from the SPATS pinger setup for depths from 80 m 
to 500 m. Both pressure and shear wave speeds were 
mapped versus depth in firn and bulk ice. This is the 
first measurement of the pressure wave speed below 
180 m depth in the bulk ice, and the first measure- 
ment of the shear wave speed in the South Pole ice. 
The resulting vertical sound speed gradient for both 
pressure and shear waves is consistent with no re- 
fraction between 250 and 500 m depth. Currently, an 
extended analysis is underway that includes longer 
distances and also a more stringent error estimate 
for the distance. 
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Fig. 4. Preliminary results for Vpressure (top) and Vghear 
(bottom) versus depth for the bulk ice down to 500 m depth. 
The best fit is indicated. 
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500 m depth region for the pressure and shear wave 
measurements. The graphs include a fit for the sound 
speed gradient. The best fit yields a radius of cur- 
vature of ^ 30 km for both the pressure and shear 
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